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There have been a number of recent spectral models that have been successful in reproducing the observed 
X-ray spectra of galactic black hole candidates (GBHC). However, there still exists controversy over such issues 
as: what are the sources of hard radiation, what is the system's geometry, is the accretion efficient or inefficient, 
etc. A potentially powerful tool for distinguishing among these possibilities, made possible by the Rossi X-ray 
Timing Explorer (RXTE), is the variability data, especially the observed phase lags and variability coherence. 
These data, in conjunction with spectral modeling, have the potential of determining physical sizes of the system, 
as well as placing strong constraints on both Compton corona and advection models. As an example, we present 
RXTE variability data of Cygnus X-l. 



1. INTRODUCTION 

Many black hole candidate (BHC) systems 
show three or more spectral states: the "low" (X- 
ray hard), "high" (X-ray soft), and "very high" 
(X-ray soft, plus a hard, power-law tail) states. 
Typically, sources are seen to be in the low state 
when they are below ~ 10% of their Eddington 
luminosity, and they tend to exist in the high and 
very high states at observed luminosities above 
10% ^Edd @- In this article, we will concentrate 
on the properties of the low state, for which we 
have obtained ~ 20 ks of Cyg X-l data with the 
Rossi X-ray Timing Explorer (RXTE) [[2|3j]. 

A variety of models have been proposed for the 
low state of BHC in general, and for the low 
state of Cyg X-l in specific. Currently, these 
models fall into two main classes: the Comp- 
ton corona models [Q~[| and the Advection Dom- 
inated Accretion Flow (ADAF) models @|l|- 
For the Compton corona models, the currently 
favored geometry is one in which a central, spher- 
ical corona is surrounded by an exterior, cold 



disk [g]|||j Hj. This is similar to the geometry 
proposed by the ADAF models [pTJ|Jll[ ; however, 
the ADAF models invoke different sources for the 
Comptonization seed photonsn, and they make 
more detailed predictions for the dynamics of the 
inner region. Specifically, the advective inner re- 
gion is moving quasi-radially, at close to the free- 
fall velocity, toward the black hole. 

We note that Cyg X-l has upon occasion, in- 
cluding very recently, transited to the high state 
113] . Different scenarios have been proposed for 
this transition in both the Compton corona pic- 
ture H and the ADAF picture Q. We will 
not discuss this transition here. Instead, we will 
concern ourselves with how variability data, i.e. 
the phase/time lags and the coherence function, 
can be used to constrain models within a given 
state. As described below, these data currently 
pose challenges for both the Compton corona and 
ADAF models of the low state. 



x ADAFs take the Comptonization seed photons to be cy- 
clo/synchrotron photons generated within the advective 
region itself. 



2. SPECTRAL MODEL OF CYG X-l 

As an example of a successful (in terms of fit- 
ting the energy spectra) low state model, we con- 
sider the spherical corona model of Dove et al. 
0,|L|l||- Its basic features are as follows. The in- 
ner region of the accretion system is modeled as a 
spherical cloud with total optical depth ^1.6 and 
an average temperature of 87 keV. This spherical 
region is then surrounded by a cold (T < 150 eV), 
geometrically thin disk which reprocesses hard ra- 
diation from the corona, as well as provides seed 
photons for Comptonization. The overall quality 
of the fit of this model to our Cyg X-l data is 
quite good (x? d = 1-6)- We will use this as our 
"straw man" model when we discuss the limits 
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that phase lags place on coronal models (§4.4) 



3. POWER SPECTRAL DENSITY 

The power spectral density (PSD) is calcu- 
lated by taking lightcurves from the data, divid- 
ing them into segments of equal length, and then 
taking the Fast Fourier Transform (FFT) of each 
data segment. The squared amplitude of each in- 
dividual FFT (for a given lightcurve) is then aver- 
aged together. Furthermore, we usually average 
over (logarithmically spaced) Fourier frequency 
bins as well. This yields the resulting PSD for 
each lightcurve. Here we choose a normalization 
such that the integral of the PSD over Fourier 
frequency yields the square of the total root mean 
square (rms) amplitude of the variability p3 . For 
a given narrow Fourier frequency interval, the rms 
is (to within a factor of v2) the fractional amount 
by which the lightcurve is sinusoidally modulated 
in that given frequency interval. 

In Figure [y we plot the PSD of the low energy 
band (0 - 3.86 keV) for our Cyg X-l data. This 
PSD, in both shape and amplitude, is very similar 
to what has been previously observed. We note, 
however, that there is a weak, broad feature at 
0.005 Hz, with rms ~ 1.5%. A Lomb-Scargle pe- 
riodogram |1J] indicates that the significance of 
this feature is ~ 50%. However, for reasons dis- 
cussed in §ra, this feature may be somewhat more 
significant than that. 
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Figure 1. Power Spectral Density of the (0-3.86 
keV) band of Cyg X-l. The solid line corresponds 
to the "effective" noise level. For the normaliza- 
tion of Miyamoto jjJI, this is 2/(Count Rate), di- 
vided by the number of FFTs and frequency bins 
averaged over. FFTs for lightcurves of 1024, 128, 
and 32 s duration were combined for this plot. 



4. PHASE LAGS 

4.1. Basics 

Fourier phase lags (and equivalently Fourier 
time lags) are calculated from the FFTs con- 
structed from data segments of two different, but 
concurrent, lightcurves. Let s(t) be a "soft en- 
ergy" light curve and h(i) be a "hard energy" 
lightcurve. The Fourier phase lag, 4>(f), is then 
just the phase of the (complex) quantity 

(S*(f)H(f)) , (1) 

where S(f) and H(f) are the FFTs of the con- 
current data segments of s(t) and h(t), respec- 
tively. The * denotes a complex conjugate and the 
angle brackets denote an average over data seg- 
ments and/or Fourier frequency bins. The time 
lag, r(/) = 0(/)/2tt/. 

In general, both the phase lag and time lag are 
non-constant functions of Fourier frequency, /. 
However, we naively expect that the longest ob- 
served time lag in the system is no longer than 
the longest distance in the system divided by the 
slowest propagation speed. As shown in Fig. @, 
this longest time for Cyg X-l is ~ 0.1 s. This is 
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Figure 2. Time lags (as a function of Fourier fre- 
quency) in Cyg X-l for the (14.09 - 100 keV) 
vs. (0 — 3.86 keV) energy bands. Open dia- 
monds, the hard lags behind the soft; * the soft 
lags behind the hard. Horizontal lines are char- 
acteristic timescales for a 10 M Q black hole. Dot- 
ted lines are the (radial) free-fall timescales from 
R = 30, 100 GM/c 2 , and solid lines are the light 
crossing times for R = 30, 100 GM/c 2 . 



considerably longer than both the expected free- 
fall timescales for ADAF models and the sound 
travel time of hot corona models. Therefore, 
these data may pose a problem for both of these 
models. 

Similarly, the shortest observed time lag should 
be longer than the smallest scale (over which sig- 
nificant luminosity is generated) divided by the 
fastest propagation speed. As shown in Fig. g, 
the shortest time lags observed in Cyg X-l are 
comparable to a light crossing time, and therefore 
may be relevant to determining physic al p arame- 
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ters for Compton corona models (cf. 

If the coherence function is unity (cf. §^j), 
one can often think of the phase lag in terms 
of a transfer function. For this case, we have 
H(f) = A(f)exp[i4>(f)]S(f), where A{f) is a 
real-valued amplitude. In the time domain, this 
means that h(t) is the convolution of s(t) with a 
linear transfer function. Specifically, 



h(t) 



t r (t — t)s(t) dr 



The Fourier transform of £ r (r) is then just 
A(/)exp[#(/)]. 

4.2. Simple Interpretations 

Imagine that we have a source of fluctuations 
that produces soft X-ray photons that is some 
distance from a region that produces hard X- 
ray photons. If the disturbances can propagate 
from the soft photon producing region to the hard 
photon producing region — without dispersion — 
then we expect there to be a time delay between 
the soft and hard photons that is independent 
of Fourier frequency, f. The time delay at all 
Fourier frequencies will simply be the distance be- 
tween the soft X-ray source and the hard X-ray 
response divided by the propagation speed. This 
means that the Fourier phase lag will increase lin- 
early with / (modulo integer multiples of 2ir). 

As shown in Fig. @, the time lags observed 
in Cyg X-l are not independent of Fourier fre- 
quency. If anything, the Fourier phase lag is more 
approximately independent of /. A crude, but il- 
lustrative, model of this is as follows. If H(f) 
is the transform of the hard photon light curve 
and S(f) is the transform of the soft photon light 
curve, let 



H(f) = Aex P (±*A0) S(f) , 



(2) 



where both A and A</> are constants, and the + is 
for / > and the — is for / < 0. [The antisymme- 
try in the relative phase is due to the fact that real 
light curves produce Fourier transforms where 
H(f) = H*(—f).) Taking the inverse transform 
of Aexp(±iA<f>), the transfer function, £ r (r), can 
be written in the time domain as 

sin(A0) 



t r {r) = A 



cos(A(p)S(i 



(3) 



That is, a fraction Acos(Acj)) of the hard variabil- 
ity is exactly coincident with the soft variability, 
while a (typically smaller) fraction is delayed from 
the soft variability with ar" 1 "tail" in the trans- 
fer function. For A ~ 1 and A</> ~ 0.1 radians, 
roughly 90% of the soft and hard lightcurves are 
exactly coincident with one another. 

4.3. Simple Propagation Models 

It is often convenient to think of propagating 
disturbances in terms of sources and responses. 
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Figure 3. Phase lags (radians) for the data pre- 
sented in Fig. y. Lines cor respond to the simple 
propagation model of § |4.3| . The solid line is for 
c p ss 0.01 c, and the dashed line is for c p « 0.1 c, 
where c is the speed of light. 



For example, we might have a disturbance, ^(x, t) 
that obeys a wave equation 



/ d 2 



\dx 



s 2 J^)*(^ f ) 



-47TpD(f,t),(4) 



where po{x,t) is the source of the disturbance. 

Imagine that the source of fluctuations, po, is 
separable in space and time. Let us also assume 
that the observed soft X-ray light curve is equal 
to the disturbance, *$!, times a weighting function, 
g s (x), integrated over the system. Likewise, let us 
assume that the observed hard X-ray light curve 
is related to \l/ via a weighting function gh(x). 

As was discussed by Vaughan & Nowak JL5| , it 
is relatively straightforward to calculate the re- 
sultant phase lags for this case. Furthermore, the 
resultant "transfer function" between the soft and 
hard lightcurves has many of the qualitative prop- 
erties discussed in §4.2. 

As a simple example, let us consider the (phe- 
nomcnological) weighting functions: 



g s (r) oc 01 6 ) * exp ( -a s — I , 

Jo J \ r . 



gh{r) oc 
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where — a s and — an are constants, and 9 repre- 
sents a step function. Here we shall take a s — 
0.31 and a,h = 0.58, so that the response of the 
soft X-rays is predominantly at r ^ 15 r$ and the 
response of the hard X-rays is predominantly at 
r Ss 8 7"o. For the case of Cyg X-l, we will take 
r = 6 GM/c 2 with M = 10 M Q . The resultant 
phase lags then depend upon c p , the propagation 
speed of the disturbances (which we take to be 
uniform, with the waves directed toward a "sink" 
at the inner disk edge) jljj]. 

As shown in Fig. pi such a model qualitatively 
reproduces the observed phase lags. However, 
to obtain quantitative agreement with the phase 
lags at low Fourier frequency (~ 0.1 Hz), a very 
slow propagation speed of c p sw 0.01 c, where c 
i s th e speed of light, is required. As discussed in 
§4.1, such a slow speed is much less than both 
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the radial velocity in ADAF models, as well as 
the sound speed for Compton corona models, and 
therefore is problematic for both models. 

4.4. Comptonization Models 

Any intrinsic time delays between soft and 
hard photons for Comptonization seed photons 
will be modified as the seed photons diffuse 
through the Compton cloud p^!8| . The greater 
the observed output energy, the longer the dif- 
fusion time through the Compton corona. As 
has been discussed by Miller |L7| and Nowak & 
Vaughan [p.8| , the minimum expected time delay 
is the difference in the diffusion time for the ob- 
served hard and observed soft photons. 

We have calculated this minimum time delay 
for the Compton corona model discussed in §|2J, 
and show this along with the data in Fig. |j. De- 
pending upon whether one identifies the ~ 10 — 30 
Hz time lags or the > 30 Hz time lags as the upper 
limit for such a time lag "shelf" , an upper limit to 
the coronal radius is seen to be R rs 45 GM/c 2 
or R w 10 GM/c 2 (M = 10 M Q ), respectively 
For either case, such a small radial extent of the 
corona poses a severe problem for ADAF mod- 
els, as shown in Fig. 0, which would then have 
great difficulty explaining the longest time lags. 
(There still remains the possibility for Compton 
corona models that the longest time lags are "in- 
trinsic" time lags from the outer cold disk, merely 
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Figure 4. Time lags (as a function of Fourier 
frequency) for Cyg X-l data. The solid lines 
correspond to the minimum expected time lag 
for a spherical Compton corona (with parame- 
ters given in §g) of radius R — 45 GM/c 2 and 
R = 10 GM/c 2 (Af = 10 M Q ). Thick line is the 
magnitude of the Poisson noise limit above 30 Hz. 



"reprocessed" by the corona |l7|]lq ].) 

To be conservative, we identify the 10 — 30 Hz 
data as the upper limit for the shelf, as the > 30 
Hz data shows a strong loss of coherence (cf. Fig 
||) . We clarify this point below where we describe 
exactly what we mean by the coherence function. 

5. COHERENCE FUNCTION 

The coherence function has been extensively 
discussed by Vaughan & Nowak p5[ . Essentially, 
it is the average, normalized amplitude of the 
cross power spectral density (which is what one 
is computing to find the time lags). Similar to 
cq. [j], the coherence function, C(f) is given by 

(S*(f)H(f)) 2 



C(f) 



(6) 



(S 2 (/))(ff 2 (/)) ' 

Various methods exist for determining, and min- 
imizing, the effects of Poisson noise on the esti- 
mate of this function pa] . 

As also discussed by Vaughan & Nowak fll5| , 
most previously considered mechanisms for pro- 
ducing variability in BHC systems (i.e. multiple 
shots or flares, any "nonlinear" processes, etc.) 
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Figure 5. Coherence function (as a function of 
Fourier frequency) for the data presented in Fig. 
H Note how remarkably close the coherence is to 
unity in the range ~ 0.02 — 10 Hz. The loss of co- 
herence at frequencies below ~ 0.02 Hz is real, as 
is the recovery to near unity coincident with the 
0.005 Hz feature from Fig. [y. The loss of coher- 
ence at high frequency (^ 30 Hz) may be influ- 
enced by uncertainties associated with the Pois- 
son noise level and the instrumental deadtime. 



will lead to a strong loss of coherence. This makes 
the near unity coherence shown in Fig. R quite re- 
markable. The near unity coherence seen between 
~ 0.02 — 10 Hz argues for a single source of dis- 
turbances and/or a global response in Cyg X-l 
over (these rather disparate) timescales pq| . 

In a very real sense, the coherence function 
is a measurement of the degree of constancy of 
the phase lags between soft and hard photons as 
one averages over individual data segments and 
Fourier frequency bins (as discussed in 



4.1). If 



the phase is constant from data segment to data 
segment, the coherence is unity. A variable phase, 
as might be caused by two or more uncorrelated 
processes relating soft and hard photons, will lead 
to coherence loss p5| . As the coherence drops 
rapidly above <; 30 Hz, we cannot be sure that 
the observed time lag in this regime is indicative 
of solely the diffusion through a Compton corona. 
For this reason, we take the 10 — 30 Hz data as 
the conservative estimate. 



Note also that there is a definite loss of co- 
herence below ~ 0.02 Hz. There is a clear trend 
which is only broken near 0.005 Hz. At that point, 
there is a recovery of coherence to near unity, 
coincident with the feature seen in the PSD (cf. 
Fig. nj). This indicates that, although the Lomb- 
Scargle periodogram indicated a significance of 
only 50%, this feature may actually be significant. 
At the very least, it argues that PSDs of Cyg X-l 
should be extended to frequencies as low as 10 -3 
Hz whenever statistics permit. 

6. SUMMARY 

We have discussed the role of Fourier 
phase/time lags and variability coherence in con- 
straining spectral models of BHC. Much of the 
theory behind this work can be found in the pa- 
pers of Dove et al. |,|, Miller 0, Nowak & 
Vaughan @, and Vaughan & Nowak @. To 
illustrate our points, we have presented RXTE 
data of Cyg X-l (cf. §|). 

Despite the general success of spectral mod- 
els, such as the Compton corona model discussed 
in §0, the variability data pose substantial chal- 
lenges. For example, both the ADAF models and 
Compton corona models (if they rely on sound 
speed propagation timescales) have substantial 
difficulty in explaining the longest observed time 
lags. Explaining the observed unity coherence in 
Cyg X-l will be difficult for all models unless they 
can postulate a global variability mechanism. 

The shortest observed time lags were also seen 
to provide upper limits to the size of a spherical 
Compton corona. The "recovery" to near unity 
coherence at ~ 0.005 Hz, coincident with a weak 
feature in the PSD, indicates the possibility of an 
extremely low frequency oscillation in Cyg X-l. 

In short, no currently proposed model for the 
low state of BHC in general, or the low state 
of Cyg X-l in specific, can claim to completely 
describe both the spectral and variability data. 
However, with the advent of RXTE and Bep- 
poSAX, which between them are capable of broad 
spectral coverage and fast timing, only models 
that can explain, or at least be consistent with, 
both the spectral and variability data should be 
considered to be viable. 
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